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Abstract. A method for calculating the short-range order part of the free 
energy of order-disorder systems is proposed. The method is based on the 
apllication of the cumulant expansion to the exact configurational entropy. 
Second-order correlation corrections to the mean-field approximation for the free 
energy are calculated for arbitrary thermodynamic phase and type of interactions. 
The resulting quadratic approximation for the correlation entropy leads to 
substantially better values of transition temperatures for the nearest-neighbour 
cubic Ising ferromagnets. 
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In a recent paper jjj we developed a method of deriving very simple polynomial 
approximations for the variational configurational entropy in which the well-known 
and reliable cluster variation method (CVM) of Kikuchi || was used as a starting 
point. It was observed that in the absolute majority of the CVM approximations 
the number of variational parameters in the entropy is substantially greater than in 
the internal energy. As a result, the minimization of the free energy with respect 
to those variables which are not present in the internal energy reduces to the 
maximization of the entropy. The corresponding equations simply relate different 
variational parameters, since the configurational entropy contains no information 
about the character and strength of interactions in the system. Most of the variables 
in the CVM free energy are therefore redundant and can in principle be eliminated 
by solving these equations. However, the equations are non-linear and apparently 
cannot be solved exactly; even if this was possible, the resulting expression for the 
entropy would be far too complicated. To overcome this difficulty, it was proposed to 
expand all the redundant parameters in powers of cumulants of those pair and higher- 
order averages which enter the expression for the internal energy; such cumulants are 
usually sufficiently small. In this case the initial equations reduce to those for the 
corresponding expansion coefficients, and the latter can be easily solved. The final 
approximation for the entropy is a polynomial of some order in the cumulants of only 
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those averages which are present in the internal energy. The proposed method can be 
regarded as a simplification of the CVM. 

The described approach, however, combines two very different approximations, 
since its aim is to find an approximate analytical solution of the already approximate 
CVM equations (which are usually treated only numerically) . In the present letter we 
show that it is possible to apply the same idea of the cumulant expansion directly to the 
exact configurational entropy, thus avoiding using the CVM for constructing the initial 
approximation. The general formalism is presented and the first non-zero (quadratic) 
corrections due to short-range order (SRO) to the mean-field approximation (MFA) for 
the free energy are calculated. Considerably more complicated calculation of higher- 
order terms is left for a discussion elsewhere g . It is shown that the obtained quadratic 
approximation for the correlation entropy significantly improves, in comparison with 
the MFA, calculated transition temperatures for the nearest-neighbour (NN) cubic 
Ising ferromagnets. 

The formally exact expression for the configurational entropy has the form 



S = -ks 



p2l...2W ]„ pH---ZiV 

M....JV lur l...N ■ 



(1) 



Here kg is the Boltzmann constant, the TV-site probability 

Pi.:n N ={pi---Pn) (2) 
is the statistical average of the product of N occupation numbers, and N is the 
number of lattice sites in the system. The occupation number p % n for a particular spin 
orientation (up or down) is equal to 1 or depending on whether the spin at the 
lattice site n has this orientation or not. It is convenient for our purposes to define the 
orientation index i as taking on values ±1. Equation (|l]) is a slightly more detailed 
version of the familiar result [Q 

S = — fcsTrplnp (3) 

where p is the density matrix and Tr denotes the trace of a matrix. 

The cumulant average j^] of a product of I variables is the Zth-order residue when 
all combinations of cumulant averages of lower order are subtracted from the actual 
average of the product. The cumulant averages of the occupation numbers are thus 
defined by the relations 

«) = «) c 

<«) = «) c «) c + («) c 

(<«) - (<) c «) c (<>c + («) c (<>c + («) c (<) 

+ («) c (<) c + ««) c 

etc., where the subscript c denotes a cumulant average. Correspondingly, the A^-site 
probability (||) is written as 

it:i x - (pi 1 ). ■ ■ ■ (p^) c + [ (p\ 1 P2 2 ) c (p'i) c ■ ■ ■ (p l £) c + • • • 

+ (p^) c ...(pT-1) (pT-\p^ 
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<K i pr> c (pM 4 > c <^ 5 > c ---<p^> c +--. 



+(pi i ... P i £) c . (5) 

Introducing notations 

P£:X = (<•••<) Ci\;X = (<...< ) c (6) 

and noticing that P^ = as given by the first of equations , we rewrite the rest 
of equations (Q) and equation (|J) as follows, 
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Finally, defining new quantities 
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we represent equations (f7t) in the form 
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where in the last equation 

Y%x—i N _ \ " rp-L—i 
^l...N — ^ni...) 

9> 2 (rii...n g ) 

Here notation {n\ . . . n q ) corresponds to summation over clusters consisting of lattice 
sites n±, . . . , n qi and the last term denotes all contributions containing products of 
three or more cumulants C^'.'.'.ri, • Insertion of the result for the TV-site probability 
(the last of equations (^)) into equation (Q) leads to the natural separation of the 
long-range order (MFA) and SRO (correlation) contributions to the configurational 
entropy 

5* = Sq -\- S\ 



S = -k B ^P^\nP^ 
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We now expand the correlation entropy S\ in powers of X % ^'"^ N . It can be shown that 
the linear term vanishes ||; this result is obvious, since the entropy is maximal in the 
absence of correlations. Substitution of equation (nil) into equation (pTd) then gives 



Si = -^ E^.-^EE E E «,^::X + °( c3 )- ( 12 ) 

»l...ijv q>2 p>2 (m...n q ) (m 1 ...rn v ) 

The next step is to use the relation 

C%:% =2-«j\...j q C ni ... nq (13) 

where C ni ,.. nq is the cumulant average of the product of the corresponding spin 
variables s„ which are linearly related to the occupation numbers and acquire values 
+ 1 for spin up and —1 for spin down, 

C ni ...n q — ( s ni ■ ■ • s ng) c ■ (14) 

The proof of equation (|l^) is not given here to save space M, but its validity can 
be directly checked, at least for low-order cumulants. Inserting equation ( p"3] ) into 
equation (||) and then the result into equation (jl^), we obtain 



IN 



5 1 = ^yy2-^+ i » y y y ■ ;: ---- / " /, ' | --- / ' ; - />i ' 

9>2p>2 (ni...n 9 > (mi...ro f ) Vi-W r m ■ ■ ■ r n, ■ ■ ■ r m f / 

x Cni...ri,C mi ... mp + O (C 3 ) . (15) 

The expression in the round brackets in equation (jl^) factorizes into the product 
of sums over orientations of individual spins. The result of each such summation 
depends on whether the corresponding lattice site belongs to both sets n±, . . . , n q and 
mi, . . . , m p , to only one of them, or to neither, and is 

3 

y.1 = (+l) + (-l) = (166) 
j 

ypi = p++p- = 1 (i6 C ) 

3 

for these three cases, respectively (in the second ca se th e probabilities in the numerator 
and denominator cancel each other). Equation ( |166| ) shows that the only surviving 
terms in equation ( [j"5| ) are those in which the sets n\, . . . , n q and m\, . . . , m p represent 
the same lattice cluster. Substituting equations ( |l6o| )-( |l6~c| ) into equation ( p^|) 
and neglecting products of three or more cumulants C ni ... n , we finally obtain the 
expression for the correlation entropy in the second-order approximation, 

S 1 = -k B y2-^+V y (P+P-...P+P-y 1 C 2 ni „, nq . (17) 

q>2 (ni...n,) 

In this order the correlation entropy is thus a linear combination of squares of the 
cumulants, i.e., contains no cross-terms. The minimization of the free energy then 
leads to vanishing of all those cumulants which are not present in the internal energy 
Equation ( |l7| ) for the correlation entropy has been obtained without any 
assumptions about the nature of a thermodynamic phase or interactions in the system. 
We now list several frequently encountered situations in which the result for the free 
energy has particularly simple form. 
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(i) In the absence of sublattices, when all lattice sites are equivalent (such as in 
the case of disordered or ferromagnetically ordered material), the point probabilities 
do not depend on the site index n, 



where m = (s n ) is the site magnetization. Equations (116) and (|T^) therefore become 



/1 + m, 1 + m 1—m, 1 — m\ 

50 = -JVAb(— In— + — In— ^— 1 (19a) 

5 1 = -ifc s ^(l-m 2 )- 9 ^ C n,..« 9 - ( 19& ) 

<?>2 (n-i_...n q ) 

(ii) If pairwise interactions are assumed in the case (i), then the only terms 
necessary to retain in equation ( [196] ) are those with q = 2. All other, higher- 
order cumulants vanish at equilibrium, since the internal energy is independent of 
corresponding averages. As a result, 

S 1 = -\k B {\-m 2 y 2 ]T C 2 nm . (20) 

(rim) 

At this stage it is convenient to use pair SRO parameters a nm instead of pair cumulants 
C nrrl ; the relation between them is ^ 

C nm = (l - to 2 ) a nm . (21) 

In terms of the SRO parameters equation (|2^) has the form 

Si = -\k B Y, <n - -\Nk B J2 Z l<*l ( 22 ) 

(nm) / 

where in the second equation Zi and a/ are the coordination number and SRO 
parameter for the coordination shell I, and the summation is performed over all 
coordination shells. 

(Hi) When only the NN spins interact in the case (ii), it is sufficient to take 
into account in the second of equations ( f22| ) only the term corresponding to the first 
coordination shell. The rest of the SRO parameters vanish at equilibrium for the same 
reason as before. The correlation entropy in this case is simply 

Si = -^Nk B Za 2 . (23) 

Here we have dropped the subscripts in the NN coordination number and in the NN 
SRO parameter. For the FCC lattice (Z = 12) this result coincides with that obtained 
in [Q within the tetrahedron CVM approximation. This means that the tetrahedron 
version of the CVM reproduces correctly the exact expansion for the configurational 
entropy up to at least second order in a. 

Finally, to test the obtained approximation for the free energy, we consider 
thermodynamics of the NN Ising ferromagnet, an example which corresponds to the 
simplest case (Hi). The Hamiltonian of the Ising model is written as 

H ^ Jnm^n^m h ^ ^n* (^^) 

(rim) n 
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Here J nm is the exchange integral equal to J > for NN sites and to zero otherwise, 
and h is the magnetic field. The internal energy, being the statistical average of the 
Hamiltonian, is therefore 

E = —NZJ [m 2 + (l -m 2 ) a] - Nhm (25) 

where the first two of e qua tions (Q) for the spin variables and equation (|2l]) were used. 
Combining equations ( [19 cj ) , (^J) an d ( p5[ ) , we get for the free energy F — E — TS (T 
is the absolute temperature) 

F = N 

V 2 

1 + m , 1 + m 1 — m , 1 — m \ 1 



— ZJm 2 — hmj — ~^ZiJ (l — m 2 ) - 



/w;T ( ——h; _____ + _____ | n _____ j + lk B TZa 2 



(26) 

The four terms in equation (|26|), left to right, are the MFA internal energy, correlation 
internal energy, MFA entropy and correlation entropy, respectively. The variational 
free energy is now a function of two minimization parameters m and a, whose 
equilibrium values are determined by the conditions 

OF dF 2 _ 

dm da 
They lead to the following equations for m and a, 

m(l — a) + n 



— tanh 



t 



Zt 



(28) 



where t = ksT/ZJ and /i = h/ZJ are the dimensionless temperature and magnetic 
field. Note that the first of equations (|_8|) differs from its well-known MFA counterpart 
by the factor 1 — a. The temperature of the second-order transition is the solution of 
the equation 

'd 2 F\ 

^) - °- (29) 

Combining this equation with the second of equations ( p8| ) for m = 0, we obtain the 
equation for t c , 

t c = l- {Zt c y X . (30) 

The ^-independent MFA result t c = 1 is obtained from this equation in the limit 
Z — > oo. Equation (|30| ) shows that the derived SRO corrections to the free energy lead 
to the decrease in t c of order _T _1 , as expected 0. Results of the numerical solution 
of equation ( ^0| ) for the three cubic lattices are given in table [lj they are compared 
with the best- known t c values ||. It is seen that the quadratic approximation for 
the correlation entropy improves the accuracy of the calculated t c from 20-30% in the 
MFA to 5-10%, and the accuracy increases with decreasing coordination number. 

To conclude, in this letter we have developed a method of calculating the 
correlation contribution to the free energy of a system which undergoes the order- 
disorder transition. The essence of this method is the expansion in powers of cumulant 
averages as applied to the exact configurational entropy. It has been shown that 
this approach leads to the construction of simple polynomial approximations for the 
entropy. First corrections to the MFA free energy (see equation ( |l7| ) for the correlation 
entropy) have been calculated without any assumptions about the thermodynamic 
phase or the character and range of the interactions in the system. Previously, the 
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Table 1. Normalized temperature t c = kgT c /ZJ of the second-order phase 
transition for the NN Ising ferromagnet on the simple cubic (SC), body- 
centered cubic (BCC), and face-centered cubic (FCC) lattices calculated using 
equation (BOI), i.e., within the quadratic approximation (E3|) for the correlation 
entropy. The best- known t c values (which can be considered as essentially exact), 
the MFA results and the corresponding ratios are shown for comparison. 



Lattice 


Z 


tc 


t MFA 
L c 


xexact 
L c 


tc / t e*act 


^MFA uexact 


SC 


6 


0.789 


1 


0.752 


1.049 


1.330 


BCC 


8 


0.854 


1 


0.794 


1.075 


1.259 


FCC 


12 


0.908 


1 


0.816 


1.113 


1.225 



same idea was utilized by the author to simplify the videly used CVM by eliminating 
redundant variables in the CVM variational entropy. However, the approach presented 
here is more consistent, since the smallness of the pair and higher-order cumulants is 
the only assumption used, while in Q this assumption was combined with the very 
different CVM approximation. The advantage of both proposed methods is that the 
entropy polynomial needs to be derived once and for all for a given Hamiltonian and 
type of order; after that thermodynamic calculations become extremely easy and can 
be performed either analytically or using a pocket calculator, at least in relatively 
simple cases. For the quadratic approximation considered here this has been shown 
by the above example of the NN Ising ferromagnet. The SRO corrections lead to 
considerably better agreement, in comparison with the MFA, of the approximate 
transition temperatures with the corresponding best-known values. Calculation of 
the higher-order terms of the cumulant expansion for the configurational entropy is 
less straightforward, though perfectly possible, and is beyond the scope of the present 
short letter; this issue will be discussed in a separate publication ||. 
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